hep-ph/0307364 
IFIC/03-40 
ZU-TH 11/03 



Neutrino properties and the decay of the lightest supersymmetric particle 



M. Hirsch 1 and W. Porod 2 

1 Astroparticle and High Energy Physics Group, IFIC - Instituto de Fisica Corpuscular, 
Edificio Institutos de Investigacidn, Apartado de Correos 22085, E-46071 Valencia, 

Espana 

2 Institut fur Theoretische Physik, Universitdt Zurich, 
CH-8057 Zurich, Switzerland 

Abstract 

Supersymmetry with broken R-parity can explain the neutrino mass squared 
differences and mixing angles observed in neutrino oscillation experiments. In the 
minimal model, where R-parity is broken only by bilinear terms, certain decay 
properties of the lightest supersymmetric particle (LSP) are correlated with neu- 
trino mixing angles. Here we consider charginos, squarks, gluinos and sneutrinos 
being the LSP and calculate their decay properties in bilinear R-parity breaking 
supersymmetry. Together with the decays of charged scalars and neutralinos cal- 
culated previously this completes the proof that bilinear R-parity breaking as the 
source of neutrino masses will be testable at future colliders. Moreover, we argue 
that in case of GMSB, the decays of the NLSP can be used to test the model. 



1 Introduction 



In supersymmetric models with R-parity violation (RpV) [1, 2, 3] the lightest supersym- 
metric particle (LSP) is unstable and decays. Thus astrophysical constraints on its nature 
[4] no longer apply and a priori any SUSY particle could be the LSP. On the other hand, 
most studies of LSP phenomenology at accelerators in RpV models up to now have con- 
centrated on a) the lightest neutralino [5, 6] or b) a charged scalar (most probably the 
right scalar tau) [7, 8] being the LSP. The purpose of the present work is to study the 
phenomenology of the remaining LSP candidates (charginos, gluinos, scalar quarks and 
scalar neutrinos) at future colliders for the case where R-parity is broken by bilinear terms 
only. 

Arguably at present the main motivation to study RpV SUSY models is the astonishing 
experimental progress in neutrino physics in the past few years. Super-K observations of 
atmospheric neutrinos [9], solar neutrino measurements by the SNO collaboration [10] 
and the reactor anti-neutrino experiment KamLAND [11] have finally established non- 
zero neutrino masses and mixings beyond any reasonable doubt. Bilinear RpV (BRpV) 
SUSY models necessarily produce Majorana neutrino masses and indeed in [12] it was 
shown that BRpV SUSY can explain current neutrino data, once 1-loop corrections are 
carefully taken into account. 

Can one test BRpV SUSY being the origin of neutrino masses and mixings? Con- 
ventional wisdoms says, if SUSY is to solve the gauge hierarchy problem, superpartners 
should be found at the next generation of colliders. In general, the only predictable 
difference between R-parity conserving and R-parity violating SUSY then would be a 
decaying LSP. Given the relatively small number of free parameters in BRpV SUSY, 
however, one can go further and from measured data on neutrino properties predict sev- 
eral decay properties of the LSP. This was shown for the case of the neutralino be- 
ing LSP in [5] and for charged scalar LSPs in [7]. In both cases, interesting relations 
between certain decay patterns and low-energy neutrino data have been found. The 
most important are: (i) The ratio BR(xi — ¥ l^oq') /BR(xi T QQ') — tan 2 ^tm and (ii) 
BRfa -> ev)/BR(f 1 -> fiv) ~ tan 2 9 sol . 

The minimal supersymmetric extension of the standard model (MSSM) [13] contains 
more than a hundred free parameters, most of which are soft SUSY breaking masses 
and phases. Supplementing the MSSM with universal mSugra boundary conditions [13] 
reduces this number to just four free parameters plus an undetermined sign (in addition 
to the standard model parameters). These are given at the grand unification (GUT) scale 
as: m , the common scalar mass, m^, the gaugino mass, and A , the common trilinear 
parameter. In addition one usually chooses tan (3 = v u / Vd and the sign of the Higgs mixing 
term as free parameters. Parameters at the electro- weak scale can then be calculated 
from RGE running, reducing considerably the available parameter space. 

In such a constrained version of the MSSM (CMSSM) one finds essentially only two 
LSP candidates, namely, the lightest neutralino and the right sleptons, in particular the 
right scalar tau if tan/5 is large. Arguably it is this theoretical prejudice why also in RpV 
models other possibilities so far have attracted little attention. 

Models which depart from strict mSugra in one way or another, however, can be found 
in the literature. Just to mention a few representative examples, there are string inspired 
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models where supersymmetry breaking is triggered not only by the dilaton fields but 
also by moduli fields [14]. In SO (10) or E(Q) models, where all neutral gauge bosons, 
except those forming Z and 7, have masses of the order of ttigut one expects additional 
D-term contributions to the sfermion mass parameters at incur [15]. This is equivalent to 
assuming non-universal values of mo for left-sleptons, right-sleptons, left-squarks and right 
squarks, giving rise e.g. to sneutrino LSPs. In gauge mediated supersymmetry breaking 
(GMSB) models [16] exists the possibility, that the gluino is the LSP [17]. In AMSB 
models [18] one can find parameter regions where the chargino is the LSP (nearly mass 
degenerate with the lightest neutralino). 

To be as general as possible, however, in this work we will not resort to any specific 
model of SUSY breaking. Instead we will simply point out in which way one has to depart 
from mSugra to obtain the corresponding LSP and then proceed to calculate its decay 
properties. To summarize our main result it can be said that independent of which SUSY 
particle is the LSP there is at least one ratio of decay branching ratios which is fixed by 
either the solar or the atmospheric neutrino angle, i.e. independent of the LSP nature, 
BRpV SUSY as the origin of neutrino masses is testable at future colliders. 

This paper is organized as follows. In the next section we will recapitulate the main 
features of the bilinear R-parity violating model. Then we will turn to the numerical 
calculations. Decays of squarks, gluinos, charginos and scalar neutrinos will be discussed 
in detail. Before concluding with a short summary, we will argue that also in case of 
GMSB [18], BRpV SUSY remains testable due to the decay patterns of the NLSP. 



2 The Model 

Bilinear R-parity breaking supersymmetry has been discussed extensively in the literature 
[1, 12]. We will therefore summarize only the main features of the model here, with 
emphasis on neutrino physics. The Lagrangian of the model is obtained by adding bilinear 
terms breaking lepton number to the MSSM superpotential: 

^BRpV = Wmssm _ £ ab£iL°;H^ , (1) 

and consistently the corresponding terms to the soft SUSY breaking potential: 

U so ft = U softi MSSM — ZabBitiLIH 1 ^ . (2) 

The latter induce vacuum expectation values Vi for the sneutrinos which are in turn 
responsible for mixing between standard model particles with supersymmetric particles: 
Higgs bosons with sleptons, charged leptons with charginos and, most importantly for the 
following considerations, neutrinos with neutralinos. 

The mixing of neutrinos with neutralinos gives rise to one massive neutrino at tree 
level. The other two neutrinos obtain masses due to loop effects [12]. Assuming that the 
heaviest neutrino obtains its mass at tree level, the main features relevant for our current 
purpose are the following. The mass of the heaviest neutrino is given by: 

(gM 1+ g>M 2 )\Af 
mu3 ~ 4Det(x°) 1 ' 

K = €iV d + IIVi. (4) 
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The atmospheric neutrino mixing angle is given by 

A 2 



tan# 



Atm 



A, 



(5) 



and the so-called CHOOZ angle by 



U, 



■\ 2 > ■ M ' 



(6) 



l 2 ~ IV 3 

The scale of the loop-induced solar mass is given by 



m b 



and the solar mixing angle by 



tan 9, 



sol 



V-i/,tree 



(7) 

(8) 
(9) 



where V u,tree is the tree level neutrino mixing matrix [12] . In the region where the condition 
(e2A2)/(e 3 A 3 ) < is fulfilled one finds that ei/e 2 ~ ei/e 2 . For a more thorough discussion 
see ref. [12]. 

In this model the neutrino spectrum is hierarchical and hence the neutrino mass scales 
coincide with the (square roots of) the mass squared differences measured in oscillation 
experiments. This implies that the R-parity violating parameters are significantly smaller 
than the R-parity conserving parameters: |ej| <C \fi and \vt\ <C v^. This feature allows 
for the possibility that all R-parity violating couplings can be expanded in terms of the 
ratios 



A,, 



or 



Ai 



A*' JDet(x°) Det (^ 



(10) 



Several examples of this kind can be found in [5, 7, 12, 19]. We have used this possibility 
for a systematic expansion of the R-parity violating couplings to obtain a semi-analytical 
understanding of the results presented in the following section. For completeness we want 
to note that in our model also gauginos and gauge bosons have R-parity violating couplings 
implying that one can clearly distinguish BRpV from a model where only trilinear R-parity 
violating couplings are present. 



3 Numerical results 

In this section we present various collider observables and their correlations with neutrino 
observables for the following LSP candidates: charginos, sneutrinos, squarks and gluino. 
The numerical results are obtained in the following way except if stated otherwise: (i) We 
create a random sample over the SUSY parameter space, using five free parameters: m , 
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m i/2, A , tan(3 and [i. Motivated by mSugra, we calculate the gaugino masses approxi- 
mately from mi/2, for the sfermion mass parameters we use mo directly at the electroweak 
scale. We have checked explicitly for several LSP sets that the latter simplification has no 
impact on our results. We then violate one of the following mSugra conditions at a time. 
Each condition is necessary but not sufficient to obtain the corresponding LSP, i.e. after 
calculating the SUSY spectrum we post-select points which have the desired candidate 
LSP. Chargino LSPs are obtained by the condition m 2 ~ (5/3) tan 2 dy/ m \-, sneutrino LSPs 
by m 2 L . <C ttIq, squark LSPs by m 2 Q ) rn 2 D) mfj <C m\ and gluino LSPs by m 3 <C m 2 ,mi. 
(ii) The R-parity violating parameters are added such, that A\ tm and A 2 ol are consistent 
with the experimental data. 

We want to stress that the correlations between low-energy and high-energy observ- 
ables shown in the following are predictions after a generous sampling over the SUSY 
parameter space. Much tighter correlations could be obtained, once information on at 
least a part of the SUSY spectrum is put in, see e.g. ref. [5] for the case of neutralinos. 
In particular information on x% Xki ^ h and H + would be important in this respect. 

3.1 Charginos 

For chargino LSPs possible final states are 

Xt -> J2 qq'vubtvi (11) 

q=d,s;q'=u,c 

Xt - n E qq;ltcc;lfbb;ltit; (12) 

q=u,d,s 

Xt - IPPk (13) 

Xt - ifVrVs (14) 

where k = e, /i, r and we sum over the three neutrino flavours as well as the lighter quark 
states u, d and s which cannot be separated experimentally. Note, that we calculate here 
the 3-body decays even if an intermediate real 2-body final state is possible by including 
the finite width of the intermediate states. These intermediate states contain in general 
a gauge boson, whose R-parity violating couplings to charginos are typically an order of 
magnitude smaller compared to the R-parity violating couplings of the virtual sfermions. 

Numerically one finds that the final state qq' J2j v j has usually the largest branching 
ratio (up to 65 %). Typical branching ratios for other final states not containing a top 
quark are in the range of several per mille to a few per-cent. Final states with top quarks 
are found to be always very small or kinematically closed, because m-+ does not exceed 
300 GeV in our numerical data sets. Moreover, the intermediate states are always off-shell 
in this case contrary to final states involving light quarks. 

As a first example how neutrino physics allows to predict observables for collider 
physics we plot in Fig. 1 r/m|+ [meV/(100 GeV) 5 ] as a function of the heaviest neutrino 
mass m U3 [eV]. We scale out a fifth power of the chargino mass to account approximately 
for the phase space of the decay. 1 

1 Although this was not discussed in [5] a similar correlation holds for the case of the lightest neutralino 
being the LSP. 
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Figure 1: Chargino decay width divided by m|+ as a function of the heaviest neutrino 
mass 



Fig. 1 shows an obvious correlation between chargino decay width and neutrino mass. 
However, there is a sizable spread in the prediction. Thus such a measurement could 
probably never compete with neutrino oscillation experiments in terms of accuracy. Nev- 
ertheless, from the current data on atmospheric neutrinos one can roughly predict, 



Eq.(15) can be considered as a consistency check for the completeness and uniqueness of 
the bilinear model as the main source of the (atmospheric) neutrino mass. Significantly 
smaller or larger widths would be a clear signal that BRpV cannot explain the neutrino 
data. Note that the width of the band gets reduced sizeably once some information on 
the SUSY spectrum is available. 

As discussed previously [5], ratios of different branching ratios of the LSP decays can 
trace information on ratios of RpV parameters. In case of the lightest chargino being the 
LSP we have found that various ratios are sensitive to ratios of Aj/Aj. Two examples 
are shown in Fig. 2. In this figure we show BR(x + — > ecc)/BR(x + —> ^cc) (to the 
left) and BR(x + —> l^cc)/BR{x + tcc) (to the right) as function of Ai/A 2 and A 2 /A 3 , 
respectively. Obviously measurements of these branching ratios would determine the 
corresponding ratio of A's to high accuracy. Somewhat worse results are obtained if one 
has to sum over the quarks of the first two generations. 

Since A2/A3 determines the atmospheric angle in BRpV, one expects the ratios dis- 
cussed above to be correlated with tan 2 6 Atm- That this is indeed the case is shown in 
Fig. 3, where we plot BR(x + — > fiqq)/BR(x + —> r <?<?) as a function of tan 2 ^Atm sum- 
ming over all quarks of the first two generations. For the currently preferred value of 
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1CT 2 1(T 1 10° 10 1 1CT 2 10 _1 10° 10 1 10 2 



(A1/A2) 2 (A2/A3) 2 

Figure 2: Ratio of branching ratios for chargino decay. To the left, BR(x + 
ecc)/BR(x + — > ij.cc) as a function of (Ai/A 2 ) 2 . To the right, BR(x + — > /icc)/BR(x 
rcc) as a function of (A 2 /A 3 ) 2 . 



Table 1: Ratio of branching ratios for chargino LSP decays as required by the consistency 
of the model. These ratios all trace the ratio A 2 /A 3 . The experimentally allowed range 
for the atmospheric neutrino angle, 0.3 < sin 2 (9 Atm) < 0.7 (at 3 o c.L), has been used to 
obtain the quoted ranges. Ratios have been sorted with respect to increasing uncertainties. 



Ratio 


lower bound 


upper bound 


Br (/ice) / Br (rcc) 


0.45 


1.4 


Br(fiqq) / Br(rgg) 


0.45 


1.4 


Br(/i2e) / Br(r2e) 


0.45 


1.4 


Br(nbb) 1 Br(rbb) 


0.46 


2.1 


Br(3/i) 1 Br(r2/i) 


0.28 


1.8 


Br(3//) / Br(3r) 


0.096 


1.4 



tan 2 6Utm = 1 one expect approximately equal branching ratios for these final states, 
nearly independent of any other parameter. 

Due to the fact that the chargino has many possible final states, one can devise various 
cross checks of the bilinear model. We have found that the following ratios are sensitive 
to Ai/A 2 only to a good approximation: BR(x + — > eqq)/BR(x + — > fxqq), BR(x + —> 
3e)/BR(x + -> /i2e) and fii?(x + -> 3e)/BR(x + -> 3/i). 

Similarly there are a number of ratios which trace very well the ratio A 2 /A 3 . Since 
the latter quantity is related to the atmospheric angle, using the currently available ex- 
perimental data on tan 2 (^ tm ), one can predict various ratios of branching ratios. Ranges 
allowed by current data are listed in Table 1. 
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Figure 3: Ratio of chargino decay branching ratios BR{\ 
tan 2 Q\tm- F° r the currently preferred value of tan 2 9jy tn , 
equal branching ratios for these final states. 



fxqq) I BR(x + —> T~qq) versus 
1 one expect approximately 



3.2 Sneutrinos 

In the case that sneutrinos are the LSPs they will decay according to 

i>i -> qq (16) 
h -> Ipl (17) 
Vi -> VjV k (18) 

Scanning over the parameter space one finds the following general features: (i) The 
main decay mode is for all sneutrinos u,i — > bb. (ii) The branching ratios for decays into 
charged leptons are O(10~ 2 ). (iii) The invisible decay mode into two neutrinos is O(10 -3 ) 
and below. Moreover, it turns out that the decay lengths of all sneutrinos have the same 
order of magnitude because (a) the largest couplings are of the form e^Y h and (b) the 
largeness of the solar mixing angle requires t\ ~ e 2 and one expects that also e 3 is of 
similar size. This implies that one has to sum over all sneutrinos if one considers the 
direct production at an e + e~ collider. Therefore we consider the following quantity: 

3 

a * BRilflJ) = ]T a(e + e- -> u s u r -> bblflj) (19) 

r,s=l 

In Fig. 4a - d we show the ratio of these observables as a function of (ei/e2) 2 and tan 2 6 so \. 
The correlations of the collider observables with neutrino physics are obvious. The range 
for the observables under study is a*BR(e/j) = O(10~ 2 )-O(l) fb, a*BR(/j,/j,) = O(0. 1) fb, 
a* BR(/j,t) = O(0.1)-O(10) fb at a 800 GeV e + e~ linear collider and unpolarized beams. 
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Larger cross sections are obtained for polarized beams with left-handed electrons and 
right-handed positrons. We want to stress again, that the spread in the correlations is 
mainly due to the unknown SUSY spectrum. It gets considerably reduced once informa- 
tion on the SUSY spectrum is plugged in. 

In principal one could tag the flavour of the sneutrino in cascade decays, e.g. in the 
decay Xi ^> l + v and then study the subsequent decay of the sneutrino into leptons. Due 
to the fact that all involved decays are two-body decays, one can distinguish the lepton 
stemming from the chargino from the ones stemming from the sneutrino by measuring 
the lepton energy except in the case where m-+ ~ 3mj>/2. In Fig. 4e and f we show 
correlations between the branching ratio of the muon sneutrino and (ei/e2) 2 and tan 2 8 so \ 
assuming that such a flavour tag can indeed be performed. As can be seen, current results 
for the solar mixing angle predict that BR(z> M — > e ± fi T ) /BR(z> M — > / u ± r T ) is in the range 
0.4 — 2 independent of the remaining SUSY parameters. 

In scenarios, where sneutrinos are the LSPs, the left charged sleptons are not much 
heavier. The difference between the masses of the charged left sleptons and sneutrinos 
is roughly given by: m~ L — m? ~ — cos2 ( 5m 2 y > 0. Depending on the mass difference 
the left sleptons decay in these scenarios either via three body decays, which conserve 
R-parity, into [20, 21] 

h -> vqq' (20) 
l L -> vvl (21) 

or via R-parity violating couplings into 

Il -> qq' (22) 
l L -> vl. (23) 

The latter decay modes give in principal rise to additional observables correlated with 
neutrino physics. However, we have found that for mass differences larger than ~ 5 GeV 
the three body decays clearly dominate. Therefore this additional information is only 
aviable if either tan/5 is small and/or if all particles have masses above > 400 GeV. 

3.3 Squarks 

Here we discuss the decays of the squarks of the first two generations as well as the decays 
of the lighter sbottom. The decays of the lighter stops have been discussed in detail in 
ref. [19] and are similar to the results for the sbottoms discussed below. 

In the case that squarks are the lightest SUSY particle, they will decay according to 

q -> qu (24) 
q - q'l (25) 

In case that the mass difference between the squarks of the first two generation is larger 
than approx. 5 GeV, the heavier ones will dominantly decay according to 

q - q'qq' (26) 
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(ei/e 2 ) 2 tan 2 O so i 



Figure 4: Various sneutrino observables as a function of (ei/e2) 2 (left column) and tan 2 9 so i 
(right column). <j* BRiljlj) is defined as Srs=i °"( e+e ~ - ► ^s^V - > bblflj). For discussion 
see text. 
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mediated mainly via gluino exchange. Here we have used the formulas given in [21]. This 
decay mode dominates the decays of di, because mg L > rriu L + 5 GeV for rriQ < 500 GeV 
and tan (3 > 3 due to the D-terms in the mass matrix. 

The important decay modes related to neutrino physics are those induced by the 
effective LiQD c R coupling which is proportional to e^. In the numerical results below we 
assume that the first two generations of squarks are mass degenerate. This assumption 
is motivated by the experimental constraints from meson physics, in particular for the 
K°-K° mixing [22]. However, we do not assume left and right squarks have the same 
mass. 

Typical examples are shown in Fig. 5, where we consider three different scenarios: 
(i) In Fig. 5a and b we show the case where ul and cl are the LSPs. One sees that 
the corresponding relations between tan 2 6^ (or (ei/e 2 ) 2 ) and ratios of branching ratios 
into charged leptons are extremely pronounced. Note that the decays into leptons clearly 
dominate, with a total branching ratio of 0.6 - 0.9 summed over all charged leptons. (ii) 
In Fig. 5c and d we show how the results are changed if one sums over the left and right 
squarks of the first two generations. As discussed above, d L and sl are not included in 
this sum because their R-parity violating decay modes are negligible. The correlations are 
somewhat worse compared to the previous case because ur (cr) decay into leptons only via 
their mixing with the corresponding left partner and the corresponding branching ratios 
are of order 10 _3 -10 -2 . The branching ratios of cLr and sr into leptons is approximately 
in the range 0.05 - 0.5. (iii) In Fig. 5e and f we show the correlations if b\ is the LSP. 
One sees that the ratio of branching ratios into top quarks is nicely related to the ratios 
of the €i squared. In the case that this decay is suppressed the decays into a lepton and a 
c-quark give similar, although somewhat worse, results. In this case the branching ratio 
is of order 10~ 2 . 

Finally, we want to note that these correlations are hardly affected by QCD corrections 
because they nearly drop out by taking the ratio of branching ratios. We have checked 
this explicitly by adopting the formulas given in [23] to the bilinear model. 



3.4 Gluinos 

In the case that the gluino is the LSP, it decays according to 

9 -> WQ (27) 

g - ^qq' (28) 

9 -> Vi9 (29) 

These decays proceed via virtual squarks. For this reason one expects also correlations 
between ratios of branching ratios into Uqq' and ratios of q and, thus, the solar mixing 
angle. We have adopted the formulas given in [24] for the calculation. The general 
features of a gluino LSP are: (i) The decay into the final state vjob dominates, where we 
sum over all neutrinos. This can be seen in Fig. 6. (ii) The sum of the branching ratios 
of the decays into l^bt (I = e, /i, r) final states is of order 10 _2 -10 _1 if summed over all 
charged leptons. Note that at the LHC O(10 5 ) gluino pairs can be produced per year if 
rrig = 500 GeV. (iii) All other decay modes are at most of order 10~ 2 . 
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Figure 5: Ratios of branching ratios for squark decays as function of (ej/ej) 2 (left column) 
and tan 2 9 soi (right column). 
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Figure 6: BR(g — > 66i/) as a function of 



The important decay modes for testing correlations between gluino branching ratios 
and neutrino angles are the final states l^bt (I = e,/x,r). The sum of these decay modes 
is normally of the order of a few per-cent. In these decays the same class of couplings 
induced by the effective vertex L{T L B C R being proportional to h b €i are probed as in the 
decays of the lighter stop [19] and in the decays of the lighter sbottom. As can be seen 
in Figs. 7a and b, there is a clear relation between the ratios of the final states into l{tb 
and the corresponding ratio of e^. In case of the ratio BR(g — > etb) / BR(g — > fitb) this 
implies a clear correlation with the solar mixing angle as can be seen in Fig. 7b. 

3.5 Gravitino 

In gauge mediated SUSY breaking [16] the gravitino G is the LSP. Depending on the scale 
of SUSY breaking, its mass is typically of the order eV up to several MeV. If R-parity is 
broken, the gravitino decays, but it is long lived from the point of view of collider physics, 
because the width is proportional to an R-parity violating coupling squared and the ratio 
{m,Q/msusY) A - This implies that the gravitino will escape the detector before decaying. 

However, one can consider the decays of the next-to- lightest SUSY particle (NLSP). 
For a gravitino with a mass of 1 eV one finds a partial width for NLSP decays into gravitino 
of O(10 -3 ) eV and significantly smaller values for larger gravitino masses. Partial widths 
of the lightest neutralino (or chargino) into R-parity violating final state are of the same 
order of magnitude. This implies that branching ratios for R-parity violating final states 
of the neutralino are at least of O(10~ 2 ). For slepton NLSPs typical partial widths into 
R-parity violating states are of 0(1) eV in case of staus and thus clearly exceed the decay 
into a gravitino. Very similar arguments apply to all the other NLSP candidates discussed 
in the previous sections. Therefore, we conclude that in GMSB models with bilinear R- 
parity violating terms, the decays of the NLSP can be used to establish the correlations 
with neutrino physics we have discussed. 
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4 Summary 

We have calculated the decay patterns of various possible LSPs within bilinear R-parity 
violating supersymmetry. The main conclusion of the present work is that whichever 
SUSY particle is the LSP, measurements of branching ratios at future accelerators will 
provide a definite test of bilinear R-parity breaking as the model of neutrino mass. In 
case of GMSB, where the gravitino is the LSP, we find that correlations with neutrino 
physics exist for the decays of the NLSP. 

One can state the above more carefully. Observation of a decaying LSP would provide 
proof that R-parity is violated. Measuring ratios of branching ratios then presents the 
ultimate cross-check of the completeness and uniqueness of the bilinear model. The most 
robust predictions of BRpV are shown in this paper, but many ratios of branching ratios 
are tightly constrained by neutrino physics. In fact, for several different LSP candidates 
many different decay channels have sizable branching ratios which should follow the spe- 
cific patterns discussed in the previous sections. Thus this simplest model of R-parity 
violation can be over constrained by the measurements we have discussed and - in this 
sense - easily ruled out. 
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